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Controller Design for Robust Output Regulation of 

Regular Linear Systems 

Lassi Paunonen* 


Abstract 

We present three dynamic error feedback controllers for robust output regulation 
of regular linear systems. These controllers are (i) a minimal order robust controller 
for exponentially stable systems (ii) an observer-based robust controller and (iii) 
a new internal model based robust controller structure. In addition, we present 
two controllers that are by construction robust with respect to predefined classes of 
perturbations. The results are illustrated with an example where we study robust 
output tracking of a sinusoidal reference signal for a two-dimensional heat equation 
with boundary control and observation. 

Keywords: Robust output regulation, regular linear systems, controller design, feed¬ 
back. 


1 Introduction 


The topic of this paper is the construction of controllers for robust output regulation of 
linear inhnite-dimensional systems. The goal in this control problem is to design a control 
law for a linear system 


x{t) = Ax{t) + Bu{t) + w{t), x(0) = Xq G X 

y{t) = Cx{t) + Du{t) 


(la) 

(lb) 


in such a way that the output y{t) converges asymptotically to a given reference signal 
yref{t) despite the external disturbance signal w{t). In addition, the controller must toler¬ 
ate small perturbations and uncertainties in the parameters {A, B, C, D) of the plant flTl) . 
The robust output regulation problem was hrst studied for hnite-dimensional systems in 
the 1970’s most notably by Francis and Wonham and Davison [^, and since then 

the theory of output regulation has been been actively developed for inhnite-dimensional 


systems j2l|, 0, [^. 22. Ill . 

The most recent developments in the held are related to the study of output reg¬ 
ulation and robust output regulation for inhnite-dimensional systems with unbounded 
input and output operators, and especially for regular linear systems jl^, 30, 2^ which 
are often encountered in the study of partial diherential equations with boundary control 
and observation [^. In particular, the characterization of the solvability of the output 
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regulation problem using the so-called regulator equations was extended for systems with 
unbounded operators B and C in E, [m |. and the internal model principle of robust 
output regulation was generalized for regular linear systems in (l^ . 

In this paper we continue the work begun in 19j. The main emphasis in the refer¬ 
ence jl^ was on studying the properties of robust controllers and on characterizing the 
solvability of the robust output regulation problem. In this paper we concentrate on 
designing actual controllers that achieve robust output regulation for the regular linear 
system m- As our main results we present three different robust controllers. Two of 
these controllers employ structures that are familiar from the control of systems with 
bounded operators B and C, and the third employs a completely new complementary 
controller structure. 

The reference signal Vrefi,-) and the disturbance signal w(-) are assumed to be generated 
by an exosystem 


v{t) = Sv{t), 
w{t) = Ev{t) 

Vrefit) = -Fv{t) 


n(0) = vq eW 


(2a) 

(2b) 

(2c) 


on a finite-dimensional space W = C^. Here S' is a matrix with eigenvalues cr(S') = 
{iui,... ,iuq} C iM. The main objective in this paper is to achieve robust output regu¬ 
lation for the system ([T]) by choosing appropriate parameters (^i, Q 21 K) for the dynamic 
error feedback controller 


z{t) = Qiz{t) ^ 2 e(f), 
u{t) = Kz{t). 


^(0) — zq Ei Z 


(3a) 

(3b) 


where e{t) = y{t) — yrefit) is the regulation error. 

The main tool in constructing robust controllers is the internal model principle, which 
provides a complete characterization of the controllers that achieve robust output regu¬ 
lation for the system ([T]) and for the reference and disturbance signals generated by the 
exosystem ([2]). In particular, this fundamental result tells us that the control problem can 
be solved by including a suitable internal model of the dynamics of the exosystem into 
the controller ([3]), and by choosing the rest of the parameters of the controller in such a 
way that the closed-loop system consisting of the plant and the controller is stable. The 
classical definition of the internal model (also referred to as the p-copy internal model) 
requires that if p is the dimension of the output space Y and if S has a Jordan block of 
dimension associated to an eigenvalue iuk, then the operator Qi must have at least p 
independent Jordan chains of length greater or equal to associated to the same eigen¬ 
value iuk [s, 16|. In this paper we also use an alternative definition for an internal model, 
called the Q-conditions, which is applicable even if Y is infinite-dimensional HQ. 

The first controller in this paper presented in Section 0] is constructed by choosing Qi 
as the internal model of the exosystem ([2]) and by stabilizing the closed-loop system with 
suitable choices of Q 2 and K. It is well-known that if the plant ([T]) is exponentially stable 
and S is diagonal, then this very simple structure is extremely effective. Indeed, this con¬ 
troller structure has been successfully used on several occasions for inhnite-dimensional 
systems with bounded and unbounded input and output operators 27, il QSQ. The 
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most important advantages of this controller structure is that due to the internal model 
principle, this controller is of minimal possible order, and if dimF < cx) then the resulting 


ipc 

controller is finite-dimensional. In [23| this type of structure was used for regular linear 
systems on Hilbert spaces and with U = Y. In this paper we present a minimal order 
controller that solves the robust output regulation for a regular linear system flTl) on a 
Banach space X, without restrictions on the input and output spaces, and with the most 
general choices for the stabilizing operators Q 2 and K. 

In Section m we in addition present a separate version of the minimal order controller 
for a situation where robustness is only required with respect to a predefined class Oq 


of perturbations. The design is motivated by the recent observation [l8|, [iSj that in 


such a situation the robust output regulation problem may be solvable with a controller 
incorporating a reduced order internal model that is strictly smaller than the full p-copy 
internal model. In particular, in [l^ such a controller was successfully designed for a 
given class Oq of perturbations. In this paper we present a new controller that solves 
the robust output regulation problem for a stable regular linear system and for a given 
class Oq- This new controller has the advantage over the one presented in 1^ in that the 


controller is of minimal order, and it is finite-dimensional whenever dimH < cx). This 
controller is new even for finite-dimensional linear systems and for infinite-dimensional 
systems with bounded operators B and C. 

The second robust controller of this paper presented in Section O employs a novel 


structure that was introduced in 15 for construction of controllers with reduced order 


internal models. In particular, the system operator Qi of the controller has a triangular 
structure that is naturally complementary to the structure of observer-based robust con¬ 
trollers j3, 10|. The main advantages of this new controller are that it has the natural 
structure for the inclusion of the p-copy internal model into the dynamics of the controller, 
and that it can robustly regulate plants that have a larger number of inputs than out¬ 
puts. The construction of this second controller is a new result even for finite-dimensional 
linear systems and for infinite-dimensional systems with bounded operators B and C. In 
Section E] we also use the same structure to generalize the original reduced order internal 


model based controller in [15(] for regular linear systems. 


Finally, the third robust controller presented in Section O is an observer-based con¬ 
troller that employs the triangular structure that was successfully used for robust output 
regulation of systems with bounded B and C together with infinite-dimensional diagonal 


exosystems in . In this paper we generalize the observer-based construction in [10| to 


regular linear systems with nondiagonal exosystems. 

As the first main result in this paper we present the internal model principle. This 
result was first generalized for regular linear systems in [l^ in the more general setting 
of infinite-dimensional exosystems and strong stability of the closed-loop system. In this 
paper we introduce it for regular linear systems with finite-dimensional exosystems and 
exponential closed-loop stability. We demonstrate that the exponential closed-loop sta¬ 
bility allows simplifying general assumptions of the theorem, and show that the regulation 
error has an exponential rate of decay. 

We illustrate the construction of controllers by considering the robust output regula¬ 
tion problem for a two-dimensional heat equation with boundary control and observation. 
We begin by stabilizing the system with negative output feedback, and we subsequently 
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construct a minimal order controller that achieves robust tracking of a sinusoidal reference 
signal. 

The paper is organized as follows. The standing assumptions on the plant, the ex¬ 
osystem and the controller are stated in Section [2J In Section [3] we formulate the robust 
output regulation problem and present the internal model principle. The minimal order 
controller for stable systems is presented in Section 01 In Sections El we introduce the 
new controller structure for robust output regulation. Finally, the observer-based robust 
controller is constructed in El The robust output tracking of the two-dimensional heat 
equation is considered in [71 


2 The Plant, The Exosystem and The Controller 


If X and Z are Banach spaces and A : X —)• Z is a linear operator, we denote by T>{A), 
Af{A) and TZiA) the domain, kernel and range of A, respectively. The space of bounded 
linear operators from X to Z is denoted by £(X, Z). If A : X —>• X, then <7(^4), crp{A) and 
p{A) denote the spectrum, the point spectrum and the resolvent set of A, respectively. 
For A G p{A) the resolvent operator is given by i?(A, A) = (A — A)~^. The inner product 
on a Hilbert space is denoted by (•, •). For an operator A : T>{A) C X —)■ Zi x • • • x Z„ 
we use the notation A = where A^ : T’(H) C X —)■ Z^, for all k G {!,... ,'n.}, 

to signify that Ax = {A^x)^^-^ for x G T>{A). On the other hand, for an operator A G 
£(Xi X • • • X X„, Z) we use the notation A = (Hi,..., H„), meaning that Ax = YJk=i AkXk 
for all X = {xk)k=i G Xi X • • • X X„. 

We consider a linear system ([T]) on a Banach space X with state x(t) G X, output 
y(t) G X, and input u(t) G U. The spaces U and Y are Hilbert spaces. The operator A : 
V{A) C X —>• X generates a strongly continuous semigroup T{t) on X. For a fixed Aq G 
p{A) we define the scale spaces Xi = (T>(H), ||(Ao — H)-||) and X_i = (X, ||i?(Ao, H)-||) 
(the completion of X with respect to the norm ||i?(Ao, H)-||) [1, Sec. H.5]. The extension 
of A to X_i is denoted by H_i : X C X_i —)■ X_i. 

Throughout the paper we assume that ([T]) is a regular linear system [2 


30, 26, 25 


In particular, B G C{U,X_i) and C G £(Xi,X) are admissible with respect to A and 
D G C{U,Y). The operator C in fllbl) is replaced with its A-extension 


Cj^x = lim XCR{X,A)x 

X^oo 


with TX{Ca) consisting of those a; G X for which the limit exists. If C G i2(X, X), then 
C\ = C. For a regular linear system we have 71{R{X, A)B) C V{Ca) for all A G p{A), 
and the transfer function of ([T]) is 26|, Sec. 4] 


P{X)=CaR{X,A)B + D VAgp(H). 


Finally, we dehne Xb = RiA) + 7l{R{Xo, A_i)B) C R(Ca), which is independent of the 
choice of Aq G p{A). 

Assumption 2.1. The pair {A,B) is exponentially stabilizable and there exists L G 
£(X, X) such that A + LCa generates an exponentially stable semigroup. 
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The stabilizability of {A, B) means that there exists K G C{Xi, U) such that {A, B, K\) 
is a regular linear system for which I is an admissible feedback operator, and {A+BK\)\x 


28 


generates an exponentially stable semigroup 

The exosystem ([2]) is a linear system on the finite-dimensional space W = for 
some r e N, and S G C{W) = E G £(1T,X), and F G CiyV^Y). We assume the 
geometric multiplicity of each of the eigenvalues ct(S') = {iukYk=i C iM. is equal to one. 
We denote by G N the size of the Jordan block associated to iuk G cr(S'). The following 
standing assumption is crucial for the solvability of the robust output regulation problem. 
An immediate consequence of this assumption is that in order to achieve robust output 
regulation it is necessary that dim 17 > dimR. 


Assumption 2.2. For every k E {I,... ,q} we have iUk G p{A) and P{iuJk) G C{U,Y) 
is surjective. 


The dynamic error feedback controller ([2]) is an abstract linear system on a Banach 
space Z. The operator Qi : ViQi) <Z Z ^ Z generates a semigroup on Z, and Q2 G 
C{Y,Z) and K G C{Zi,U) is admissible with respect to Qi. The operator A in ([3]) is 
replaced with its A-extension K\. 

The closed-loop system consisting of the plant ([T]) and the controller ([3]) on the Banach 
space Xe = X X Z with state Xeit) = z{t)Y' is of the form 

Xe{t) = A^X^t) + BeV{t), Xe(0) = XeO G Xg, 
e{t) = CeAXe{t) DeV{t), 


where e{t) = y{t) — yrefif) is the regulation error, Xeo = (xo,^o)^, C'e = (C'a, DKa), 

( A_i BKa \ 


De = F 


B. = 


' E ^ 


[G2Ca Gi+G2DKa) 

The operator Ag : T>{Ag) (Z Xg ^ Xg has the domain 

V{Ag) = { (x, z)eXbX V{GGj I A.^x + BKaz G X } 

where Xb = V{A) + A(R(Ao, A_i)R), and V{Cg) = V{Ca) x V{Ka) D V{Ag), Bg G 
C{W,X X Z) and Dg G C{W,Y). Here CgA is the A-extension of Cg. 


Theorem 2.3. The closed-loop system [Ag, Bg, Cg, Dg) is a regular linear system. 


Proof. See jl9|. Sec. 8]. 


□ 


3 The Robust Output Regulation Problem and the 
Internal Model Principle 

We can now formulate the robust output regulation problem. We consider perturbations 
{A, B, C, D, E, F) E O of {A, B, C, D, E, F) where the operators in the class O of admis¬ 
sible perturbations are such that (i) the perturbed plant {A, B, C, D) is a regular linear 
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system and (ii) iuk G p{A) for all k G {1,... ,g}. These two conditions are in particu¬ 
lar satished for all bounded and sufficiently small perturbations to {A, B,C, D), and for 
arbitrary bounded perturbations to the operators E and F. 

The Robust Output Regulation Problem. Choose the controller (^ 1 ,^ 2 ,-R) in such 
a way that the following are satisfied: 

(a) The closed-loop semigroup T^if) is exponentially stable. 


(b) For all initial states Xeo G and Vq & W the regulation error satishes e“'e(-) G 
L^(0, cx); Y) for some a > 0. 


(c) If the operators {A, B, C, D, E, F) are perturbed to {A, B, C, D, E,F) G O in such 
a way that the closed-loop system remains exponentially stable, then for all initial 
states Xeo G and Vq eW the regulation error satishes e“'e(-) G oo;Y) for 
some d > 0. 


We have from 1^, Sec. 4] that for initial states Xgo G E>{Ae) the regulation error 
e(-) is a continuous function and hmi_,.oo e(f) = 0 whenever the property (b) holds. Thus 
for such initial states the condition e“'e(-) G L‘^{0,oo-,Y) for an a > 0 implies that the 
regulation error decays to zero at an exponential rate. 

In the following we present two dehnitions for an internal model (l^ . 19l |. In Dehni- 
tion l3. II “independent Jordan chains” refer to chains originating from linearly independent 
eigenvectors of Qi. 


Definition 3.1. Assume dimR < cx). A controller (^ 1 ,^ 2 , A') is said to incorporate a 
p-copy internal model of the exosystem S' if for all A; G {1,..., g} we have 


dim.N’iioJk — Qi) > dim Y 


and Qi has at least dimR independent Jordan chains of length greater than or equal 
to Uk associated to the eigenvalue iojk- 

Definition 3.2. A controller {Qi,Q 2 ,X) is said to satisfy the Q-conditions if 

TZ{iuJk - Gi) ^ 71 ( 02 ) = { 0 } Vfc G {1,... ,g}, (4a) 

M{52) = {0}, (4b) 

C - C/i) VA: e { 1 ,...(4c) 

The following lemma gives a sufficient condition for invariance of the ^-conditions in 
the situation where the matrix S of the exosystem is diagonal. 

Lemma 3.3. Let S be a diagonal matrix. If the operators (^ 1 ,^ 2 ) satisfy the G-condi- 
tions, and if K : V{Gi) <Z Z ^ Y is such that Af{iuk — Gi) C M{K) for a// /c G {1,..., q\, 
then also {Gi + G 2 K,G 2 ) satisfy the G-conditions. 

Proof. Since S is a diagonal matrix, we have = 1 for all k G {!,..., g} and the 
condition (Hell is trivially satisfied. Because the condition N'{G2) = { 0 } is identical for 
both (^1 -1-^2^', G2) and (Gi, G2), it is sufficient to show that IZ{iuJk — Gi+G2K) r\Il(G2) = 
{ 0 } for all k. To this end, let w = (iuk — Gi — G2K)z = G2y for some k G { 1 , ..., g}, 
z G BiGi) and y E Y. This implies {iUk — Gi)z = ^2(2/ + Xz) E IZ{iUk — Gi) H Il{G2), 
and we thus have z G Mifuk — Gi)- Due to our assumptions we then also have Kz = 0 
and w = {iujk — Gi)z = G2y, which hnally imply tc = 0 due to (Ha)) . □ 
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The following theorem presents the internal model principle for regular linear systems 
with hnite-dimensional exosystems and exponential closed-loop stability. 

Theorem 3.4. Assume that the controller stabilizes the closed-loop system exponentially. 
Then the controller solves the robust output regulation problem if and only if it satisfies 
the Q-conditions. 

Moreover, if dimY < oo, then the controller solves the robust output regulation prob¬ 
lem if and only if it incorporates a p-copy internal model of the exosystem. 


Proof. Since generates an exponentially stable semigroup and S' is a matrix with 
spectrum on zR, the Sylvester equation SS = A^T^ -\- has a unique solution S G 
C{W,Xe) satisfying 77(S) C T>{Aff) [^. Because an exponentially stable semigroup is 
also strongly stable, and since zR C p(^e), we have from [l9|, Thm. 7.2] that the controller 
satisfies the ^-conditions if and only if it solves the robust output regulation problem as 


defined in the reference [l^. The definition of the robust output regulation problem 


m 


can be obtained from our problem statement with the following modifications: 

(i) The exponential closed-loop stability is replaced by strong stability. 

(ii) It is assumed that for all admissible perturbations the Sylvester equation ES' = 
AeT + Be has a solution. 


(hi) The condition e" e(-) G T^(0, cx); Y) for Xeo G Xe is replaced by hmi_j.oo e(-) = 0 for 
XeO G 'D{Ae). 

We begin by showing that under the assumption of exponential closed-loop stability 
the two conditions in (iii) are equivalent. We prove this only for the nominal closed- 
loop system {Ae, Be,Ce, De). For perturbed parameters the situation can be handled 
analogously. We have from [l^, Lem. 4.3] that 


Xeit) = Te{t)XeO - Te{t)Yvo + TTsit)vo 


for all Xeo e Xe and Vq G W, and since {Ae, Be, Ce, De) is a regular linear system, 


e{t) = CeATe{t){XeO - Yvq) + (CeS De)Ts{t)vo 


is defined for almost all f > 0. In addition, if Xeo G B{Ae), then e{t) is continuous 
and is given by the above formula for all f > 0. The error contains the two terms 
e{t) = ei{t) -\- e2{t). The second term e^i-) = (CgE -|- De)Ts{-)vo is continuous and it 
is either nonvanishing or identically zero jl9l . Lem. A.l]. Since Te{t) is exponentially 
stable, for some a > 0 the hrst term satisfies e“'ei(-) G oo;Y) for all Xeo G Xe 

and Vo G W. These properties imply, under the assumption of exponential stability of 
the closed-loop system, that the regulation error satisfies e“e(-) G oc;Y) for all 

Xeo G Xe and vq & W ii and only if hmi_,.oo e(f) = 0 for all Xeo G V(Ae) and vq G W. 
Thus the conditions in (iii) are equivalent. 

Assume now that the controller satisfies the ^-conditions. The class of admissible 
perturbations in this paper is strictly smaller than the class of perturbations in 


19 


because exponential stability is stronger than strong stability, and because TS = AeT+Be 
has a solution for any perturbations for which the closed-loop system is exponentially 
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stable [20|. Because of this, and because we assumed the exponential closed-loop stability, 
we have from [l^, Thm. 7.2] that the controller satisfying the ^-conditions solves the 
robust output regulation problem as dehned in this paper. 

Conversely, we can now assume that the controller solves the robust output regulation 
problem. It then follows from the proof of 1^, Thm. 7.2] that the controller must satisfy 
the ^-conditions provided that the class of admissible perturbations contains = 0 
(corresponding to the zero disturbance signal) and arbitrary bounded perturbations to 
the operator F of the exosystem. Because these perturbations do not affect the stability 
of the closed-loop system, they also belong to the class O of perturbations in this paper. 
This concludes that the controller indeed satisfies the ^-conditions. 

Finally, if dimR < oo, we similarly have from (l^ . Thm. 6.2] that the controller solves 
the robust output regulation problem if and only if it incorporates a p-copy internal model 
of the exosystem. □ 


4 The Minimal Robust Controller for Stable Systems 

In this section we construct a minimal order robust controller under the assumption that 
the system operator A of the regular linear system ([T]) generates an exponentially stable 
semigroup and the matrix S of the exosystem is diagonal, i.e., 

S = diag(ia;i, ia; 2 ,..., iojq) G 

We begin by choosing the parameters of the controller. In this controller structure the 
system operator Qi contains precisely the internal model of the exosystem ([2]). This is 
achieved by defining Z = and 

Gi = dmg{tujY ,e £(Z), K = eKo = e{Kl ...,Kl)e C{Z, U) 

where £ > 0 and Kq G jC{Z, U). We choose the components Kq G jC(Y, U) of Kq in such 
a way that the operators P{iUk)KQ are invertible. This is possible due to the assumption 
of surjectivity of P{iuk), and can be achieved, for example, by choosing Kq = P^iukY 
(the Moore-Penrose pseudoinverse of P{iuk)) for all k G {1, ... ,q}. Finally, we choose 

^2 = {G^2)1=i = i-iPGu:k)K^r)l=i e CiY,Z). 

If we make the choice Kq = P^iUkY, then G 2 = ~W for all /c G {1,..., q}. 

Theorem 4.1. Assume that the semigroup T{t) generated by A is exponentially stable 
and S is a diagonal matrix. Then there exists £* > 0 such that for any 0 < e < e* 
the controller with the above choices of parameters solves the robust output regulation 
problem. 

In particular, the operators {Gi, G 2 ) satisfy the G-conditions and the closed-loop system 
is exponentially stable for all 0 < e < e*. 

Proof. We begin by showing that the controller satisfies the ^-conditions. Since Kq were 
chosen in such a way that P{iuk)KQ are invertible for all A; G {1,..., q}, we have that 
Af{G 2 ) = {0}. Let fc G {1,..., g}, zi G Z, and y E Y he such that 2 = {iuk — Gi)zi = 
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Q 2 y- The diagonal structure of Qi implies that we then necessarily have 0 = — 

— {P{iuk)KQ)*y, which is only possible if y = 0 since P{iuk)KQ is invertible. This further 
implies = ^ 2 ?/ = 0. Since k E {I,... ,q} and 2 ; G Tl{iuJk — Qi) H N'{Q 2 ) were arbitrary, 
this concludes Tl{ioJk — Gi) H TZ{G 2 ) = {0}. Finally, since = 1 for all /c G {1,..., q}, 
the condition flTc|l is trivially satisfied. 

We define H = {Hi, H 2 ,..., Hq) G C{Z, X) by choosing 

Hk = R{tu;k,A_i)BK^ 

for all k G {1,... ,g}. Due to the diagonal structure of Gi, it is easy to see that this 
operator is the unique solution of the Sylvester equation HGi = A_iH + BKq. Clearly 
TZiH) C Xb and we can dehne Cq = CaH + DKq G C{Z, Y). The operator Cq is of the 
form Co = (Cq, ..., Cq). A direct computation shows that 

Co" = CaHj, + DK^ = CARitUk, A_i)BK^ + DK^ = P{iUk)K^, 

and thus Cq = — ^ 2 - 

It remains to show that there exists e* > 0 such that the closed-loop system is 
exponentially stable for all 0 < e < £*. The closed-loop system operator is given by 

S, fel°DK) ■ = ( <"• X ^ I Xl-i- + e X } , 

If we choose a similarity transformation 

Q. = C/ "'f) = Q-y e C{x X z) 


we can define A^ = QeAf,Q~^ with domain V{Ae) = {xg G Xg | Q~^Xe G V{Ae) }. Using 
R{H) C Xb and 7l{A_iH + BKo) = TZ{HGi) C X the condition Q~^Xe G V{Af,) for 
Xe = (x, z) E X X Z becomes 


Qg ^Xg G V(Ae) ^ 


j —X -|- eHz E Xb 
^ — A—ix eA—iHz Y eBKqz E X 

Xg G D(x 4 ) X Z, 


and thus D(Ag) = P{A) x Z. Now for any Xg = (x, z) E P{Ag) a direct computation 
using HGi = A_iH + BKq and CaH + DKq = —G2 shows that 


x4.gXg f^gxdg 


—X 


eHz^ 


({A- eHG2Ca)x + e{-A_iH - BKq + HGi)z - e^HG2G2Z 

y —G2CAX + {Gi — eG 2 G^z 

[{A - eHG 2 Ca)x - e^HG 2 G* 2 z\ 
y —G2CAX + {Gi — £^2^2)^ / 

'(A-EHG2CA 0 \,2(^-HG2G^\ 

\ -G2CA Gi-eG2G*2) lyO 0 ) 
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Since C is admissible with respect to A, we have from the results in [1, Sec. III.S.c] that 
there exists > 0 such that A + eHQ 2 C\ generates an exponentially stable semigroup 
provided that 0 < e < £i. Moreover, Lemma lA.ll shows that the semigroup generated 
by Qi — £^ 2^2 is exponentially stable for all £ > 0, since 

are invertible for all k G Because Ca is an admissible input operator for 

A — eHQ 2 CA, Q 2 G £(y,Z), and the diagonal operators generate exponentially stable 
semigroups, the semigroup generated by the triangular operator is exponentially stable 
for all 0 < £ < £ 1 . Furthermore, because the second term is a bounded operator, it 
follows from standard perturbation theory of semigroups and similarity that there there 
exists £* > 0 such that A^. is exponentially stable for all 0 < £ < £*. 

Since the controller satisfies the ^-conditions and the closed-loop system is exponen¬ 
tially stable for all 0 < e < £*, we have from Theorem Id.41 that for any 0 < £ < £* the 
controller solves the robust output regulation problem. □ 

Remark 4.2. The controller presented in this section can also be used if the plant 
is initially unstable but can be stabilized with output feedback, i.e., there exists an 
admissible feedback element Ki G C{Y, U) such that the semigroup generated by {A + 
BKi{I — DKi)~^Ca)\x is exponentially stable. Indeed, in such a case the controller can 
be designed for the stabilized system ((A -|- BKi{I — DKi)~^Ca)\x^ B{I — KiD)~^, (/ — 
DKi)~^C\, (/ — DKi)~^D). This procedure is demonstrated in Section [3 

Remark 4.3. If the plant is real in the sense that P{—iu) = P{iu) for all a; G M, if 
Y = U = C”^, and if the exosystem is of the form 

S = diag(ia;i, —iui ,..., iuq, —iuq, 0) G 

then (^ 1 , Q 21 K) can be chosen to be real matrices. Indeed, in this case we can choose 

ft = diag (g!, ..., G?, Oj.p) , Gf = 

K = e{K ^,..., where = {ReP{iuJk)\ImPiiukY) G for A; G {1,..., q} 

and = P(0)1 G and dually Q 2 = (^ 2 )^ ^ where G 2 = ^ M^pxp 

k G {1,... ,g} and = —ly G The controller incorporates a p-copy internal 

model of the exosystem, and if we apply a unitary similarity transformation 

Q = diag(Qo, • • •, Qo, dy), ^ ^ [ily -i/y) ’ 

then {Q*QiQ, Q*G 2 , KQ) coincides with the controller constructed in the beginning of this 
section. From this it follows that there exists e* > 0 such that the closed-loop system is 
exponentially stable and the real controller solves the robust output regulation problem 
for all 0 < £ < e*. 

4.1 Controller With a Reduced Order Internal Model 

In this section we construct a minimal order controller for a version of the robust output 
regulation problem where the controller is only required to tolerate uncertainties from a 
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given class Oq of admissible perturbations 0,0. More precisely, in part (c) of the robust 
output regulation problem we only consider perturbations such that (A, B, C, D, E, F) G 
Oo and for which the perturbed closed-loop system is exponentially stable. We again 
assume that the plant is exponentially stable, the matrix S is diagonal, and we in addition 
assume that P{ioJk) are boundedly invertible for all /c G {1,..., q}. 

The class Oq in the control problem can be chosen freely, but it is assumed that all 
perturbations (^1, R, C, £), E, F) in Oq are such that (i) the perturbed plant (A, B, C, D) 
is a regular linear system and (ii) iuk G p{A) and the transfer function P{iuk) = 
CAR{iUk, A)B + D is boundedly invertible for all k G g}. Both of these re¬ 

quirements are in particular satisfied for sufficiently small bounded perturbations of A, 
B, C, and D. Being given such a class Oq, we begin the construction of the controller by 
defining 


Sk — 


span|P(ia;fc) ^{CR{iuk, A)Eek + Fck) {A, B, C, D, E, F) G C>o} 


C U 


for k G {1,... ,g}, where {ek)k=i is the Euclidean basis oi W = O. We further dehne 
Pk = dimiSfc. The controller that we construct contains a reduced order internal model 
where the number of copies of each of the frequencies iuk of the exosystem is exactly pk- 
It should be noted that this controller differs from the minimial order controller with a 
full internal model only in the situation where pk < dimR for at least one k G {1,..., q). 

Define Z = Yi x ■ ■ ■ x Yg where W = if pk < dim Y and Yk = Y if pk = dim Y or 
Pk = oo. We choose 

Gi = diag^zwi Jy,,..., toogly,) e C{Z), K = eKo = e(Kl..., K^) G £(Z, U) 


where £ > 0 and Kq G C(Yk, U) are such that 


^0 


{ul,.. if Pk < dimF 

P{iuk)~^ if Pk = dimE or pk = oo 


where {Uk}^!^ C t/ is a basis of the subspace Sk- Finally, we choose 

g2 = {-{P{tUk)K^or)Li^^Y,Z). 

For those k E {1,... ,q} for which pk = dimF or pk = oo we then have Q 2 — ~Iy- 

Theorem 4.4. Assume that the semigroup T(t) generated by A is exponentially stable, 
S = diag(za;i,... , iug), and P{iujk) are boundedly invertible for all k E {1,... ,q}. Then 
there exists £* > 0 such that for any 0 < e < £* the controller with the above choices of 
parameters solves the robust output regulation problem for the class Oq of perturbations. 

Proof. Let £ > 0, {A, B, C, D, E, F) E Oq and k G {1,..., g} be arbitrary and denote 
Pk = —CR{iuk, A)Eek — Fck. We begin by showing that there exists z G M{iujk — Gi) 
such that P{iuk)Kz = pk. If k is such that Yk = Y, we can choose z = (zi, ..., Zg) E Z 
such that zi = 0 for I ^ k and Zk = ^P{iuJk)P{iuJk)~^yk- Then clearly z G Mifuk — Gi) 
and P{iUk)Kz = eP{iUk)KQZk = P{iuJk)P{ioJk)~^yk = Uk- R remains to consider the 
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situation pk < dimF. Since P{iujk) ^Vk £ Sk by definition, and since {ul,... is a 
basis of Sk, there exist C C such that 

Pk 

P{iuk)~^yk = 

1=1 

Choose 2 ; = [zi,... ,Zq) such that zi = Q ioi I ^ k and Zk = ^ Then 

clearly 2 ; G Af{iuJk — Gi) and 


Pk 


P{iuk)Kz = eP{iuk)KQZk = P{iuk)^aiu[ = P{iuk)P{iuk) ^yk = yk- 


1=1 


Since k G {1 ,... ,q} and {A, B, C, D, E, F) G Oq were arbitrary, we have from 19|, Thm. 


5.1] that the controller solves the robust output regulation problem for the class Oq of 
perturbations if the closed-loop system is exponentially stable (see the proof of Theo¬ 
rem [23D- 

It remains to show that there exists £* > 0 such that for every 0 < £ < £* the closed- 
loop system is exponentially stable. However, if we define H = (ifi,..., Hq) G X) 
by choosing Hk = R{iuk, A_i)BKq, then H is the solution of the Sylvester equation 
HGi = A_iH + BKq, and the stability closed-loop system can be established exactly as 
in the proof of Theorem 14.11 since we again have C\H + DKq = —Q^- □ 


5 The New Robust Controller Structure 


In this section we introduce the new controller structure for robust output regulation of 
linear systems. This controller has the natural structure for the inclusion of a p-copy 
internal model into the dynamics of the controller. The construction of the controller is 
completed in steps. Some of the choices of the parameters require certain properties from 
the associated operators, and these properties are verified in Theorem 15.11 We begin by 
assuming that dim F < 00 . The case of an infinite-dimensional output space is considered 
separately for a diagonal exosystem in Section 15.11 

Step 1°: We begin by choosing the state space of the controller as Z = Zq x X, and 
choosing the general structure of the operators (Gi, G 2 , K) as 

. _ (G, G2 {Cj, + DK^) \ _ (G2 \ 

0 + BK^ + L{Ca + DK^))' \L j ' 


and K = (i^i, —K 2 ). The operator Gi is the internal model of the exosystem ([2]), and 
it is defined by choosing Zq = F"^i x • • • x F""^, and 

Gi = diag(jf,..., G C{Zq), K, = {kI,..., K1) g £(Zo, U). 

Here for each G {1,..., g} we have 


/ ii^kly Iy 


\ 


Jr 


Y 


iUkW 


Iy 

iuJklY/ 


G £(F"'=), 


(5) 


V 
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and K’I = ..., G C{Y^’‘, U), where G N is the dimension of the Jordan 

block in S associated to the eigenvalue iuk G cr^S). We choose the components G 
C(Y, U) of each in such a way that P{iuk)K^^ G C{Y) are boundedly invertible. This 
is possible since P{iuk) are surjective by assumption, and can be achieved, for example, 
by choosing = PliukY for all k G {1,... ,g}. For Z > 2 we can choose freely, 
e.g., = 0. 

Step 2°: By Assumption 12.11 we can choose K 2 G £(Xi,Z7) and Li G £(F, X) in such 
a way that (A_i + (here is the A-extension of K 2 ) and A + LiC\ generate 

exponentially stable semigroups. For A G p{A + LiC\) we define 


Pl{X) = CaR{\, A_i + L,Ca){B + L,D) + D. 


The identity PiiioJk) = {I — C\R{iUk, A)Li)~^P{iuk) and the choice of Ki imply that 
PL{ioJk)Ki^ G C{Y) are boundedly invertible for all k G {1,... ,g}. The domain of the 
operator Qi is chosen as 


TX{Gi) — I G Zq X Xb 


A.iXi + BK^Xi G X } . 


Step 3°: We define H = (iJi, FTa, • • •, H^) G £(Zo, X) where Hk = {Hi Hi ..., H^^l G 
C{Y^^,X) and 

Hi = j2{-ly-^R{^ook,A_, + L,CAy+^-lB + L,D)Kl. 

i=i 

We have from js^. Sec. 7] that {A + LiCa, B + LiD, D) is a regular linear system, 
and the resolvent identity in js^. Prop. 6.6] implies RiH) C Xb- We can therefore define 
Ci = CaH + DK, eC{Zo,Y). 

Step 4°: We choose G 2 G C{Y, Zq) in such a way that the semigroup generated by 
G 1 + G 2 C 1 G C{Zo) is exponentially stable (i.e., the matrix is Hurwitz). The detectability 
of the pair {Gi, Gi) is proved in Theorem 15.11 below. Finally, we define L = Li + HG 2 G 
C{Y,Z). 

Theorem 5.1. Assume dimX < 00 . The controller with the above choices of parameters 
solves the robust output regulation problem. 

In particular, the controller {Qi,Q 2 ,K) has the following properties: 

(i) The operator Qi generates a semigroup on Z and the controller {Qi,Q 2 ,K) incor¬ 
porates a p-copy internal model of the exosystem. 

(ii) The operator H is the unigue solution of the Sylvester eguation 

HG, = {A_, + L^Ga)H+ {B + L,D)K^, (6) 

and the pair (C'i,Gi) where Gi = GaH + DKi G £(Zo,Y) is exponentially de¬ 
tectable. 

(hi) The semigroup generated by Af. is exponentially stable. 
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Proof. We begin by proving part (i). We have that 


01 = 


0 



where ( ^ ) and ( g ^ ) are admissible with respect to ^ ^. It is now straightforward 

to use the results in [30, Sec. 7] to verify that 0i with the proposed domain generates a 


strongly continuous semigroup on Z = Zq x X. Moreover, it is easy to show that K\ = 
K = (Ri, —K^). For every A; G {1,..., g} the matrix Gi clearly satisfies dim.N'iiuk — 
Gi) = dim Y = p and it has exactly p Jordan blocks of size x associated to iuk- Due 
to the triangular structure of Oi, the controller therefore incorporates a p-copy internal 
model of the exosystem. 

We will now show that H is the solution of the Sylvester equation ([2]). Denote 
Al = A_i + LiCa and = B + LiD for brevity. Due to the structure of the operator 
Gi it is straigtforward to see that an operator H G £(Zo,X) such that 'TGyH) C VifJfC) 
is the solution of HGi = A^H + B^Ki if and only if for all k G {1,..., g} we have 

{lUk - Al)HI = BlK’I^ 

Ouk - Al)HI + Hi = BlK’I^ 


{lUk - Al)HI’^ + 77 ^”' = 

where H = (i7i,... , Hq), and Hk = (77^,..., 77^'“) G X). For each 7 G {1,..., g} 

the above system of equations has a unique solution 

Hi = jz{-l)'HR{iUk,AL)^+^-^BLK\P 

i=i 


Thus 77 defined in Step 3° is the unique solution of ([S]). 

We will now show that (Gi, Gi) is exponentially detectable. We can do this by showing 
that for all 7 G g} and 2 ; G M{iuJk — Gi) with 2 ; 7 ^ 0 we have Giz 7 ^ 0 [l2|, Thm. 

6.2-5]. To this end, let 7 G {1,..., g} and 2 ; G Af{iuk — Gi) such that 2 ; 7 ^ 0 be arbitrary. 
From the structure of Gi we have that 2 ; = {zi,... ,Zq) where 2 ;^ = 0 for Z 7 ^ 7, and further 
Zk = (zl, 0,..., 0) G R"'''. Using Hi = R{iujk, AlJBlK’I^ we see that 

Gi;2 = Ga77^ + DKiz = Gj,HkZk + DK^zu = G^Hlzl + 7J77f 4 = PL{tUk)K^hl + 0 


since zl 7 ^ 0 , and since we chose 77^ iii such a way that P{iUk)K'l^ and PL{iUk)K^^ are 
boundedly invertible. 

It remains to show that the closed-loop system is exponentially stable. With the 
chosen controller (^ 1 ,^ 2 , 77) the operator A^ becomes 


Tie 


/ 7l_i BKi -BK^ 

GsGa Gi + G2D771 G2GA 

V LGa LDKi 7 l_i + BK^ + LG a 
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with domain T>[Ae) equal to 








A-iX + BKiZi — BK 2 X 1 G X 
A_iXi + BK2X1 E X 


If we choose a similarity transform E C{X x Zq x X) 


Qe 


(I 0 
0 I 
\-I H 




we can define = Q^A^Q^ ^ on X x Zq x X. If we denote Xg = {x, Zi, Xi) E X x Zq x X 
and use TZi^H) <Z Xb, the domain of the operator Ag satishes 


i>(A) = { 


— \ Xe ^ X X Zc) X X 


Qg ^Xg E 


= \^Xg E Xb X ZqX Xb Qg ^Xg 


V{A ,)} 
e V{A ,)}. 


For Xg = (x, zi^xi) E Xb x Zq x Xb we thus have 


Q;^Xg E V{Ag) ^ 


A_ix + BKiZi — BK^i—x + Hzi — Xi) E X 
{A_i + BK^)[-x + Hzi-xi) E X 

(A_i + BK^)x + B{Ki - K^H)zi + BK^x^ E X 
BKiZ\ T A—\Hz\ — A—iXi E X 

(A.i + BK^)x + B{Ki - K^H)zi + BK^xi E X 
xi E V{A) 


since equation ([ 2 ]) implies A_iHzi + BKiz = HGiZi — Li{C^H + DKi)zi E X. The 
above conditions also imply x E Xb, and thus 


V{Ag) = jxe E X X Zq X V^A) 


(A_1 + BK^)x + B{Ki - K^H)zi + BK^Xi E x}. 


For Xg = (x, Zi, Xi) E T>{Ag) a direct computation using L = L 1 +G 2 H, Gi = GaH+DKi, 
and HGiZi = (^-1 + TiC'A)Ff^i + {B + LiD)KiZi yields 



QgAg 


X 

Zi 


\—x + Hzi — Xi^ 


f{A_i + BK^)x + B{Ki - K^H)zi + BK^xi 
I (Gi + G2{CaH + DKi)^zi — G2G\Xi 
\ (^-1 + LiGiAxi 



-1 + BK^ 

B{Ki - K^H) 

BK^ \ 

( X 


0 

Gi -\- G 2 G 1 

-G._Ca 

Z\ 

V 

0 

0 

A -|- LiGaJ 



The operator G 2 E C{Y, Zq) was chosen in such a way that G 1 + G 2 G 1 E G(Zq) is Hurwitz. 
Since (A_i ABK^jjx and A + LiGa generate exponentially stable semigroups, since B is 
an admissible input operator for (A + BK^jjx, Ga and are admissible input operators 
for A + LiGa, and Ki — K^H and G 2 are bounded, we have that the semigroup generated 
by Ag is exponentially stable, and due to similarity, the same is true for Ag. We thus 
conclude that the closed-loop system is exponentially stable. 

Because the controller incorporates a p-copy internal model of the exosystem and the 
closed-loop system is exponentially stable, we have from Theorem 15.41 that the controller 
solves the robust output regulation problem. □ 
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5.1 Controller for a Diagonal Exosystem 

In this section we consider the situation where the output space Y is allowed to be 
inhnite-dimensional and the matrix S in the exosystem is diagonal. We will show that in 
this situation the robust output regulation problem can be solved with particularly simple 
choice for the parameter G 2 of the controller. For a diagonal matrix S = diag(ia;i,..., iug) 
we choose Zq = and the internal model (Gi, Ki) of the exosystem is defined as 

Gi = diag(za;i/y,..., iugly), Ki = (iC|,..., Kf) 

where are chosen in such a way that P{iUk)K^ are boundedly invertible for all k G 
{1,..., g}. The following is the main result of this section. 

Theorem 5.2. Assume S = diag(ia;i,..., iuq). If the other parameters of the controller 
are chosen as in the beginning of Section O and if we choose 

G 2 = (G^)Li = e /:(F,Zo), 

then the controller solves the robust output regulation problem. 

If we choose = Piiiujky = P{iuky{I — GAR{iuJk, A)Li) for all k, then G 2 = —ly 
for all k. 

Proof. Since Ufc = 1 for all fc G {1,..., q}, we have H = (ifi,..., Hg) G C{Zo, X), where 
Hk = R{iujk, A-i + LiGa){B + LiD)K^. Because of this, the operator Gi = (Gj,..., Gf) 
satisfies 

Gf = GaH^ + DKl = PL{tuJk)Kl 

for all A; G {1,..., g}, which shows that G2 = —G^. The last claim of the theorem follows 
immediately from PiiiUk) = (/ — GARifaJk, A)Li)~^P{iuk)- The same identity and the 
fact that iCf were chosen so that P{iuk)Kf are boundedly invertible imply that the 
components Gf of G2 are boundedly invertible for all A: G g}. We thus have from 

Lemma lA.ll that the semigroup generated by Gi + G2G1 = Gi — G 2 G 2 is exponentially 
stable. The exponential stability of the closed-loop system can now be shown exactly as 
in the proof of Theorem 16.11 

Due to the fact that Y may be infinite-dimensional, we cannot use the concept of p- 
copy internal model. Instead, we will verify that the controller satishes the ^-conditions. 
For this we will in particular use Lemma 13.31 

Since S is diagonal, the condition fHcD is trivially satisfied. The components Gf = 
— {PL{ioJk)K\)* of G2 = {G 2 )k=i are boundedly invertible for all k G {l,...,g}. This 
implies 7V’(G2) = {0}, and also further shows that A/'(^2) = {0}. Moreover, if for some 
A: G {1,..., g} the elements (z, x) G Z, {w, v) G P(Gi) with w = (wk)l=i G Zq, and y eY 
are such that 

/iook — Gi 0 \ / w\ _ (G2\ 

\ 0 iUk-{A_r + BK^))\v)~\L)y^ 

then we in particular have z = {iuk — Gi)w = G2?/ and G^y = {ioJk — iuk)wk = 0. The 
invertibility of Gf implies y = 0 and {z, x) = ^2?/ = 0. Since A; G {1,..., g} was arbitrary. 
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this shows that the operators 



0 


, ^ 2 ) satisfy the ^-conditions. Since 


fG, G2 (Ca + BK^) 

[0 A_i + BK^ + L{Ca + DK^) 



0 

+ BK^ 



where for any /c G {1,..., g} we have Af{iujk ~ (*0^ A-i+bk^ )) C Zq x {0} C A/’((0, Ca + 
DK 2 )), Lemma [23] shows that the operators (^1,^2) satisfy the ^-conditions as well. 

Since the controller satisfies the ^-conditions and the closed-loop system is exponen¬ 
tially stable, we have from Theorem Id.41 that the controller solves the robust output 
regulation problem. □ 


5.2 Controller with a Reduced Order Internal Model 


It was shown in [15(] that the triangular structure used in this section is ideal for controllers 
with reduced order internal models. Indeed, if the internal model {Gi,Ki) is replaced 
with an appropriate reduced order internal model, the controller will solve the robust 
output regulation problem for a given class Oq of perturbations. As the final result in 


this section we present a generalization of the controller introduced in [15| for regular 


linear systems with diagonal exosystems. For this purpose we again assume that P{iuk) 
are invertible for all fc G {1,..., q}. 

Let Oq be a class of admissible perturbations. Similarly as in Section 14.11 we define 
Zq = Yi X ■ ■ ■ X Yg, and 


Gi = diag(iuJilYi ,..., zugly^) G C{Zo), K, = {kI..., Kl) G £(Zo, U) 
where G CiY^, U) are such that 


K 


k 

1 


{ul ,..., if pk < dim Y 
P{iuk)~^ if Pk = dimF oi pk = co 


in the notation of Section ITTl Moreover, we define G 2 = (—(PL(ia;fc)iFi G £(F, Z). 

The rest of the parameters of the controller (Qi, Q 2 , A') are chosen as in the beginning of 
Section 0 

Theorem 5.3. Assume S = diag{iui,... ,i^q) o,nd P{iuk) are invertible for all k G 
{!,..., g}. Then the controller with the above choices of parameters solves the robust 
output regulation problem for the class Oq of perturbations. 

Proof. If (A, B, C, D, E, F) G Oq and k G {1,..., g}, and if we choose z as in the proof 
of Theorem 14.41 (for e = 1), then it is easy to see that P{iuJk)K = i/k and G 
MfioJk — Qi). By jl9l . Thm. 5.1] the controller solves the robust output regulation problem 
for the class Oq of perturbations provided that the closed-loop system is exponentially 
stable. The stability of the closed-loop system can be shown exactly as in the proof of 
Theorem 15.21 □ 












6. The Observer-Based Robust Controller 


18 


6 The Observer-Based Robust Controller 


The observer-based robust controller structure presented in this section is based on the 
controller Hamalainen and Pohjolainen [l^ for systems with bounded input and output 
operators. The construction of the controller is again completed in steps and its properties 
are given in Theorem Ib.ll For this controller structure it is necessary to assume that the 
spaces U and Y are isomorphic. We begin by assuming that the plant has the same hnite 
number of inputs and outputs, that is, U = Y = C^. The case of an inhnite-dimensional 
output space is again considered separately for a diagonal exosystem in Theorem 16.21 

Step 1°: We begin by choosing the state space of the controller a,s Z = Zq x X, and 
choosing 

G = ( ^ 

\{B + LD)K^ 2l_i + BK 2 + L{Ci, + DK 2 )) ’ 

Q 2 = K = {Ki, K^)- The operators {Gi,G 2 ) make up the internal model of 

the exosystem ([2]) , and they are defined by choosing Zq = Y'^^ x • • • x R, and 

G, = diag(jy.... Jl') e C{Z„), G2 = (G|)t, e C{Y, Z„). 

Here are as in ([SD and G 2 = G C{Y, F”'') for all /c G {1,..., q}, where G N 

is the dimension of the Jordan block in S associated to the eigenvalue iuk G cr{S). We 
choose the components G^*’ G CiY) of each G 2 to be boundedly invertible (e.g., it is 
possible to choose G^*’ = ly for every k G 

Step 2°: By Assumption 12.11 we can choose K 21 G C{Xi,U) and L G C{Y,X) in such 
a way that (A_i -|- BK 2 i)\x (here K 21 is the A-extension of K 21 ) and A + LG a generate 
exponentially stable semigroups. For A G p{A_i + BK 21 ) we define 

Pk{\) = (C'a + DK^^)R{\, A_i + BK^^)B + D. 


Since P{iuJk) were assumed to be surjective for all k G {1,... ,g} and since U = Y = 
O, the identity PxiioJk) = P{i^k){.I — A_i)B)~^ implies that RxiioJk) are 

boundedly invertible for all A; G {1,..., g}. 

Step 3°: We define an operator H : Vi^H) C W_i —)• Zq in such a way that H = {Hk)k=i 
and Hk = where 


rik 

Hi = Y.i-iy~'Gl^iG^+DK^^)R{iujk,A_, + BK^^y+^-K 
3=1 

Since we have from js^, Sec. 7] that {A+BK21, B, Gf^+DK^i, D) is a regular linear system 
and Xb C V{G\) fl V{K2i), it is immediate that H G C{X, Zq) and 77(i?) C V{H), and 
we can thus dehne Bi = HB -|- G2D G C{U, Zq). 

Step 4°: We choose the operator Ki G C{Zq, U) in such a way that the semigroup 
generated by Gi -|- BiKi G C{Zq) is exponentially stable (i.e., the matrix is Hurwitz). 
The stabilizability of the pair (Gi, Bi) is shown in Theorem 16.11 below. Finally, we define 
K2 = K21 + KiH ^ P{.X, U) and choose the domain of the operator Qi as 


P{Gi) — G Fq X Xb 


A.ixi + B{KiZi + K^xi) G X}. 
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Theorem 6.1. Assume U = Y = C^. The controller with the above choices of parameters 
solves the robust output regulation problem. 

In particular, the controller (Qi, Q 2 , K) has the following properties: 

(i) The operator QI generates a semigroup on Z and the controller {Qi,Q 2 ,K) satisfies 
the Q-conditions in Definition ld.2L 

(ii) The operator H is the unigue solution of the Sylvester equation 

GiH = H{A_^ + BK^^) + G2 {Ca + DK^f) (7) 

on V{Ca) n V{K 2 i). Moreover, {Gi, Bi) where Bi = HB -\- G 2 D E C{U,Zq) is 
exponentially stabilizable. 

(hi) The semigroup generated by A^, is exponentially stable. 

Proof. The property that Qi with the given domain generates a strongly continuous 
semigroup can be seen analogously as in the proof of Theorem 15.11 

We will now show that H defined in Step 3° is the solution of ([7]). Denote A^ = 
A_i BK 21 and Gk = Ga + DK 21 for brevity. The structure of Gi implies that an 
operator H is the solution of GiH = HAx + G 2 GK if and only li H = with 

Hk = for all k, and for all fc G {1,..., g} we have 

- Ak) + Hi = GfGK 


- Ak) + 77”'“ = 

- Ak) = G’T’^Gk 

on V{Ga) PH{K 2 i). For every k E {I,... ,q} the above system of equations has a unique 
solution which is exactly 77^ in step 3°. 

We will now show that the pair {Gi,Bi) with 77i = HB G 2 D is exponentially 

stabilizable. This is equivalent to the pair {Bl,Gl) being exponentially detectable. Let 

k E {1,..., g} and 2; = (^i ,..., Zg) E Af{—iuk — G^). Then zi = 0 for I ^ k, and 

= (0,..., 0, z'l’^) with z^^ E Y. For any u E U we have 

{u,Blz) = {B,u,z) = ((77”''7? + Gf''77)u,4'=) = {G’^^>‘{GKR{toOk, Ak)B + D)u, z^’^) 

which immediately implies that we can have B^z = 0 only if zl'° = 0 due to the fact that 
and PKiioJk) are invertible. Since this also implies z = f) and since 7 G {1,..., g} was 
arbitrary, we have that the pair (Gi,77i) is exponentially stabilizable 0 Thm. 6.2-5]. 
Because of this it is possible to choose Ki in such a way that Gi + BiKi is Hurwitz. 

We will now show that the closed-loop system is exponentially stable. When the 
controller (^1,^2, 77) is chosen as suggested, we have that 

/ 7l_i BKi BKf: \ 

Ae = I G 2 GA Gi -|- G 2 DK 1 G 2 DK 2 1 

\-LGa BKi 7l_i + BK^ + LG a) 
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with domain 


V{A,) 


(x, Zi, Xi) G Xb ^ Zq X Xi 


A-iX + BKiZi + BK2X1 G X 1 

BKiZi + (R_i + BK^)xi G X j ■ 


If we choose a similarity transform Qe G £(X x Zq x X) 


(-1 0 0 \ 

Qe=\H I 0 = Q-\ 

0 l) 


we can dehne Ae = QeAeQ^ ^ on X x Zq x X. If we denote Xe = {x,zi,xi) E X x ZqX X, 
we have 

T^{Ae) = I Xg G X X Zq X X Qe G 'D(Ae) | 

= I Xg G Xb X Zq X Xb Qg ^Xg G 'D(Ag) | 

and for Xg = {x, zi^xi) G Xb x Zq x Xb we have 


Qg ^Xg G V{Ag) ^ 


—A_ix + BKi{Hx + Zi) + BK^i—x -|- Xi) G X 
BKiIhx + zi) + (R_i + BK^){-x + xi) G X 

(R_i + 5X3^1 )a; - BKiZi - BK^xi G X 
xi G T>(R) 


where we have used = X^ — KiH. Since the above condition also implies x E Xb, 
the domain of Ag becomes 

V{Ag) = [xgEXxZoX V{A) I (R_i + BK^^)x - BKiZi - BK^Xi E x}. 


For any Xg = {x,Zi,Xi) E T>{Ag) a direct computation using X^ = + KiH, Bi = 

HB + G 2 D, and GiHx = X(R_i + BKA)x + G2 (C'a + T>X 2^i)a; yields 

( -X \ ( {A_i +BKA)x-BK^Z i-BK^Xi 

igXg = QgAg Xx + = (Gi + {HB + G 2 DKi))zi + {HB + G 2 T>)X 2 ^a;i 

\-x + xi) \ (R_i -h LGa)xi 

A_i + BKA -BKi -BK^ \ 

0 Gi + XiXi BiK^ Li 

0 0 A—I LG a) \Xi) 



The operator Xi G C{Zq, U) was chosen in such a way that Gi-l-XiXi G C{Zq) is Hurwitz. 
Since {A_i + BKA)\x and A + LG a generate exponentially stable semigroups, since B 
and K 2 are admissible with respect to {A + XX^)|x and (R_i -|- LCa)|x, respectively, 
and since Xi and Bi are bounded, the semigroup generated by Ag is exponentially stable, 
and because of similarity, the same is also true for Ag. This concludes that the closed-loop 
system is exponentially stable. 

It remains to show that the controller satishes the ^-conditions. We begin by showing 
that {Gi,G 2 ) satisfy the ^-conditions. We have J\f{G 2 ) = {0} since G^'^ are boundedly 
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invertible for all A; G {1,, q}. If ^ G TZ{ioJk — Gi) fl Ti{G 2 ) for some k G {1,..., g}, 
there exist zi and y such that 2 ; = {iuk — Gi)zi = G 2 I/. Here 2 ; = [zi,...,Zq) with 
Zk = {z\,... ,z^'‘) G R”'', and structure of Gi implies that necessarily z'l'^ = 0. On the 
other hand, we have 0 = z^'° = which implies g = 0 since G^”'' is invertible, and 

thus z = G 2 I/ = 0. This concludes that Tl{ioJk — Gi) fl Ti{G 2 ) = {0}. Finally, a direct 
computation can be used to verify that MiioJk — GiY^~^ = {z = (zi ,..., Zg) \ = 

0, z; = 0 for / 7 ^ fc } C Ti{iuk — Gi). This concludes that (Gi,G 2 ) satisfy the ^-condi¬ 
tions. Moreover, the surjectivity of the operators G^’’ implies Zq = TZ{iuk — Gi)+7 l{G2 ), 
and we thus have Zq = TZ{iUk — Gi) © Ti{G 2 )- 

We will now show that (^ 1 ,^ 2 ) satisfy the ^-conditions. The condition ^[{^ 2 ) = {0} 
follows immediately from 7 V’(G 2 ) = {0}. If {z^x) G TZ{iuJk — Gi) O 7 ^(^ 2 ), there exist 
{zi^xi) G and y eY such that 

/iujk — Gi 0 \ \ 

BlK, A_, + BlK^ + LgJ[xJ ~[-L)y- 

where we have denoted Bi = B + LD. The hrst line implies ^ G Af{iujk — Gi) fl TZ{G 2 ) = 
{0}, and since Af{G 2 ) = {0}, we have y = 0. Thus {z,x) = ^ 2 ?/ = 0 and we conclude 
that 7l{iujk — Gi) n TZ{G 2 ) = { 0 }. 

Finally, let {z^x) G J\f{iuk — Gi)"'^~^- Since the closed-loop is exponentially stable, 
we have Z = TZ{iuJk — Gi) + Tl{G 2 ) for all A: G {1,..., q} [l^, Lem. 5.7]. Thus there exist 
{zi^xi) G T>{Gi) and y eY such that 

— Gi 0 ^ I \ 

BlK, A_i + BlK^ + LGa) 

We will show that y = 0, which will conclude that {z,x) E TZ{iuJk — Gi)- From the above 
equation we see that = {iuk — Gi)zi + G 2 y- The property {z,x) G Af{iuJk — Gi)"'^~^ 
and the triangular structure of Gi imply 2 ; G MGoJk — GiY’^~^ C TZ{iuJk — Gi). However, 
since Zq = 7l{iuk — Gi) © 77 (G 2 ), in the decomposition ^ = {iuk — Gij^^i + G 2 I/ we must 
then necessarily have G 2 y = 0, which further implies y = 0 due to J\f{G 2 ) = {0}. Since 
{z,x) E Af{iujk — was arbitrary, we have that flTc|l is satished. 

Since the controller satisfies the ^-conditions and the closed-loop system is exponen¬ 
tially stable, we have from Theorem 15.41 that the controller solves the robust output 
regulation problem. □ 

Finally, we consider the situation where Y is infinite-dimensional and the matrix S in 
the exosystem is diagonal. We choose Zq = Y^ and the internal model in the controller 
is of the form 

Gi = diag(zn;iJy,..., lUgly) E C{Zq), G 2 = (G^)Li e £(F, Zq) 

where the components G 2 are chosen to be boundedly invertible for all A: G {1,..., g}. 

Theorem 6.2. Assume S = diag(za;i,..., icUg) and P{iujk) E C{U,Y) are boundedly 
invertible for all k G {1,..., g}. If the other parameters of the controller are chosen as 
in the beginning of Section E] and if we choose 

K, = {-{G\PkGuji))\ ..., -{GlPKGujg))*) E C{Zq, U), 
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then the controller solves the robust output regulation problem. 

IfG2 = ((/ - then K, = {-ly, ..., - Jy). 

Proof. To show that the controller solves the robust output regulation problem, it is 
sufficient to show that the closed-loop system is exponentially stable, because the property 
that the controller satishes the ^-conditions and all the other properties considered in 
the proof of Theorem 16.11 remain valid for a general Hilbert space Y. 

Since Gi = diag(ia;fc/y)^^]^, the operator Bi is of the form Bi = G C{U,Y^), 

where for all A; G {1,..., q} we have 

B’^ = HkB + G’^D = G’^PK{zuk). 

since = G^iG/i + DK 2 i)P{iu!k, A_i + BK^i). This shows that Ki = —B^. The 
last claim of the theorem follows from PK{ioJk) = P{i(^k){I — K^iRiiook, A_i)B)~^, and 
the invertibility of G 2 and P{iuk) imply that Bf are boundedly invertible for all k G 
{!,..., q}. We thus have from Lemma lA.ll that the operator Gi + BiKi = Gi — BiBl 
generates an exponentially stable semigroup. The exponential stability of the closed-loop 
system can now be shown as in the proof of Theorem 16.11 □ 

7 Robust Control of a 2D Heat Equation 

In this section we consider robust output regulation for a two-dimensional heat equation 
with boundary control and observation. Set-point regulation without the robustness 
requirement was considered for the same system in B Ex. VI.2]. 

We study the heat equation 

x(^,0) = xo(0 

on the unit square f = (.^ 1 , ^ 2 ) £ = [0,1] x [0,1]. The control and observation are located 

on the parts Ti and r 2 of the boundary dVt^ where Ti = {.^ = (^i,0) | 0 < .^1 < 1/2} 
and r 2 = {^ = (^ 1 ,1) I 1/2 < < 1}. We denote Tq = dVt \ (Ti U r 2 ). The boundary 

control and the additional boundary conditions are dehned as 

nx nx nx 

^(e,A)|r,=«i(t), ^{fP)\v. = U2{tf ^(e,t)lro = 0 

for u(t) = {ui{t)^U 2 (t)) E U = C^. The outputs y{t) = (|/i(f),j/ 2 (f)) G H of the 

system are dehned as averages of the value of x{f,t) over the parts Ti and r 2 of the 
boundary, i.e., 

/•1/2 rl 

yi{t) = 2 x{fi,0;t)dfi, y 2 {t) = 2 x{fi,l-,t)dfi. 

Jo Jl/2 

We dehne Aq = A with domain V^Aq) = { x G | |^ = 0 on dfl }. We have from [1, 

Cor. 1] that with the above control and observation, the heat equation is a regular linear 
system (Hq) -B, C, D) with D = 0. The system becomes exponentially stable with negative 
output feedback, u = —kGx -|- u where k > 0 (cf. B, Ex. VI.2]). We choose k = 1, and 
dehne A = ((Ho)-i — BG)\x- 
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Our aim is to design a minimal order controller for the stabilized system {A, B, C, 0) to 
achieve robust output tracking of the reference signal y refit) = (—1, cos(7rf)). To this end, 
we choose the exosystem as IE = C^, S = diag(—ivT, 0, in), E = 0, and “ ( 1/2 ”o^ 1 / 2 ) ■ 
The reference signal Vrefit) is then generated with the choice vq = (1,1,1) of the initial 
state of the exosystem. 

Since p = dimE = 2, the internal model and the parameters of the minimal order 
controller are given by 

= diag(—ivr, —in, 0, 0, in, in) E K = eKq = e (^K^, Kq, E 

where £ > 0 and are to be chosen in such a way that the matrices P{—in)KQ, 
P{0 )Kq and P{in)KQ are nonsingular. We choose = Pi—in)~^, Kq = P{0)~^, and 
Kq = Piin)~^. Finally, Q 2 = (^ 2)^=1 where Q 2 = “-^ 2 x 2 for k E {!,..., 3}. We have 
from Theorem 14.11 that for small values of £ > 0 the controller achieves asymptotic 
tracking of the reference signal prefi'), and the control structure is robust with respect 
to perturbations in (A, B,C,0) that preserve the property {0, iivr} C piA)) and the 
exponential stability of the closed-loop system. In particular, this includes small bounded 
perturbations to the operators A, B, C, and D = 0. 

The robust controller also tolerates small perturbations and inaccuracies in the pa¬ 
rameters K and G 2 of the controller (although robustness with respect to these operators 
is not required in the statement of the robust output regulation problem). Because of this 
property, we can use approximations for the values Pi±in)~^ and P(0)“^ in Kq. In this 
example we use a truncated eigenfunction expansion of Aq in approximating the matrices 
F(0) and P(±i7r). Finally, the parameter £ > 0 needs to be chosen in such a way that 
the closed-loop is stable. 

The solution of the controlled heat equation can be approximated numerically using 
the truncated eigenfunction expansion of the operator Aq. For the simulation, the pa¬ 
rameter £ is chosen to be £ = 1/4. Figure [T] depicts the simulated behaviour of the two 
outputs of the plant. The solution of the controlled partial differential equation at time 
f = 16 is plotted in Figure El 





1 


Figure 1: Outputs yii-) and ?/ 2 (-) of the Figure 2: State of the controlled system 
controlled system. at time t = 16. 
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A Appendix 

Lemma A.l. Let Gi = diag(ia;i/y,... yiuqly) € C{Y^) and G 2 = (^ 2)^=1 E C{U,Y'^) 
where U and Y are Hilbert spaces. If the components G 2 of G 2 are boundedly invertible 
for all k E {1,..., q}, then the semigroup generated by Gi —G 2 G 2 is exponentially stable. 

Proof. Since G 1 —G 2 G 2 is a bounded operator, it is sufficient to show that ct(G'i—G 2 G 2 ) C 
C“. Since Gi generates a contraction semigroup, the same is true for Gi — G 2 G 2 , and 
thus a{Gi — G 2 G 2 ) C C“. It therefore remains to show that zM C p(Gi — (^ 2 ^ 2 )- 

Let iu E zM be such that u ^ Uk for all fc e {1,..., q}. We then have iu E p{Gi). 
If / + G^Riito, Gi)G 2 is boundedly invertible, then the Woodbury formula implies that 
iu E p{Gi + G 2 G 2 ). However, this is true since G 2 R{iu,Gi)G 2 is bounded and skew- 
adjoint. 

It remains to consider the case where iu = iun for some n E {1,... ,g}. We will 
show ||(za;„ — Gi -|- ^ 2 ^ 2 )^!! > c\\z\\ for some constant c > 0 and for all z E Z. If this 
is not true, there exists a sequence {zk)keN C Z such that ||zfc|| = 1 for all fc G N and 
\\{iun — Gi + G 2 Gf)zk\\ —)• 0 as A; —)■ 00 . Since iUn — Gi is skew-adjoint, we have 

II {iOJn — Gi -|- G2G2)^fc II > I {{iUn — Gi -|- G2G2)^fc; | 

> |Re((za;„ — Gi -|- G 2 G 2 )^fc, ^fc)) | = IIG^^fclP, 

and thus HG^^fcH —)■ 0 as fc —)■ cxo. For every k E N denote Zk = zl + zl where z\ E 
IZ{iUn — Gi), zl E J\f{iUn — Gi), and 1 = \\zk\f = Ik^lP 4- || 2 :fc|p. There exists Ci > 0 such 
that ||(za;„ — Gi)^^|| > Ci|| 2 ;^|| for all k eN. Thus 

Cl 114II < ll(*^^n - Gi)4ll = ||(*w„ - Gi)2;fc|| 

^ IK^I^n — Gi -|- G2G2)2:fc|| + IIG 2 II ||G2^fc|| —)■ 0 

as A; —>• 00 . Moreover, IIG 24 II > II(^2)~^II~^II4II) ^'Cid 

ll(<^4 4 4411 < 11^*2411 < ll<^2^fcll +11^2411 0 

as A: —)■ cx). We have now shown that zl ^ 0 and 4 0) but this contradicts the 

assumption that ^Zk\f = 1 for all k eN, and thus the original claim holds. In particular 
iUn ^ crp(Gi -|- G 2 G 2 ) and the range of iUn — Gi -|- G 2 G 2 is closed. Finally, the Mean 
Ergodic Theorem [l|. Sec. 4.3] implies that the range of iUn — Gi -|- G 2 G 2 is dense, and 
iUn E p(Gi -|- G 2 G 2 ). HI 
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